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Abstract 

The flavour singlet axial Ward identity with Osterwalder-Seiler twisted mass Wil- 
son fermion action is studied on a finite lattice, with finite fermion mass and the Wilson 
parameter r up to 1. Approach to the infinite volume chiral limit and emergence of 
the anomaly is significantly better than that obtained with ff{a) and ff{a^) improved 
fermion actions. We have shown explicitly that up to ^{g^), parity violating terms 
cancel in the Ward identity even at finite volume and finite lattice spacing. 
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1. Introduction 

It is well known that the naive Wilson fermion (V\ reproduces the chiral anomaly 
in lattice QCD in the continuum limit which was demonstrated analytically in weak 
coupling lattice perturbation theory by Karsten and Smit [2] who also showed that the 
result is independent of the Wilson parameter r. Subsequently, Kerler |3 1 starting from 
the flavor singlet axial Ward Identity also showed the emergence of the anomaly term 
in the continuum limit for small values of r. 

The Wilson fermions have a dimension five chiral symmetry breaking term which 
leads to complications in lattice simulations, such as additive quark mass renormal- 
ization, axial current renormalization and nontrivial mixing between operators. Even 
though the €^{a) lattice artifacts associated with Wilson fermions go to zero in the con- 
tinuum limit, the technical difficulties are significant in actual simulations which are 
performed in a finite volume with lattice spacing and quark mass both being non-zero. 

Currently, twisted mass Wilson QCD (tm WQCD) (for reviews, see, for example, 
131 ID IH) has become a popular way of simulating lattice QCD. In addition to the 
absence of unphysical zero modes of the Wilson-Dkac operator, the ability to overcome 
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some of the lattice renormalization problems and the automatic 0(ci) improvement are 
its obvious advantages. 

In the early days of lattice QCD, Osterwalder and Seller Q introduced a chirally 
twisted Wilson term (OStm WQCD) to avoid doublers. Compared to OStm WQCD, tm 
WQCD has a nontrivial flavour dependence which leads to flavour and parity breaking 
in the lattice theory which is expected to disappear in the continuum limit. On the 
other hand, parity violating effects of OStm WQCD may survive the continuum limit 
non-pert urbatively. 

In this work we explicitly show that the parity violating terms are absent up to 
G{g^) in the flavour singlet axial Ward Identity on the finite lattice. We study the 
dependence of the chiral anomaly with OStm WQCD on the Wilson parameter r as a 
function of the lattice quark mass am and show the independence as one approaches 
the chiral limit. We also study the volume dependence of the result as one approaches 
the chiral limit. 

The chiral anomaly on the lattice with finite volume and non-zero quark mass ac- 
tually provides an excellent laboratory to compare the OStm WQCD with unimproved 
as well as i^(fl) and ^(a^) improved actions with Wilson fermions. Our detailed in- 
vestigation provides a quantitative measure of the effectiveness of the twisted mass 
Wilson fermions in further reducing finite lattice spacing artifacts compared to &{a) 
and 0'{a^) Il8]|9][l3 improved Wilson fermions. 

Two comments, however, are in order: (i) The weak coupling perturbative analysis 
is not sensitive to some lattice artifacts that may be present in numerical simulations 
(for example, see fTTI), (ii) the positivity property of the quark determinant is not satis- 
fied for the OStm WQCD action 1 12| and hence numerical simulation with dynamical 
quarks with this action is not feasible. Numerical simulation is nevertheless possible 
with a mixed action (the OStm WQCD action for the valence quarks and the tm WQCD 
action for the sea quarks). 



2. Osterwalder-Seiler Chirally Twisted Wilson Term 

The standard Wilson fermion action 

x,y x,y 

where 
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with the variations ^ v^:^ = [1 - iysOCx] Wx, Vx^ w'x^ Vx [1 ^ '75 "x] lead to the 
flavor singlet axial Ward Identity 



(A|^75^(x)) = 2m{\i/^Ys\i/x) + iXx) 



(3) 
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where {G) denotes the functional average of & . Explanation of other terms are as 
follows: 

1 



The backward derivative, A*/(x) = - \f (x) — f [x — 



and (Xx) 



Trace[75 (GW+WG)] . 



(4) 



The Green's function 



G(x,y) 



1 



\y) 



Following the method of Kerler [3 1, to i^(g^), one arrives at 

(Zv> = 2/ e^vp;i ^VvW^>aW ^^^L'^°^(/'M«)cos(/?vfl)cos(ppfl) 
X Wo(;')[cos(/?;tfl)[m + Wb(/?)]- \r's:\x?{pxd)\{%{p)f , 



(5) 



e^vpA Ffiv{x)Fpx{x)I{am,r,L) 



(6) 
(7) 



Here 



^o(p) = -l^[(l-cos(flp^))], 



(8) 



ExpUcitly, — ( j)'*Ln,,n2,n3,n4 where ni ,n2,n3,n4 = 0, 1,2, 3, • • • . In the infinite vol- 
ume chiral limit, 7^-1. In all our plots it is the function I{am,r,L) which we have 
plotted. 

Next we consider the chirally twisted Wilson term introduced by Osterwalder and 
Seller Q which amounts to the replacement W ^ R = —iysW |I3|. (Note that Seller and 
Stamatescu |T3| introduced the generalization —iys exp(— /Sys) so that the chiral 
angle 9 can be directly related to the theta vacua of QCD. This generalization is not 
considered in this work.) 

The OStm WQCD leads to the Flavor Singlet Axial Wai'd Identity 



^ -2m Trace 75 G'" ~ {x'" (x)) 

= -2m Trace 75G"" - Trace 75 {G'"R + RG"') 



with 
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where 



withVi = -a^^[D^,D^] and Vf 



(11) 



Evaluation of {x'^ix)) following a similar method as above now produces some parity 
violating terms which are not of the form FF and therefore do not contribute to the 
anomaly, however, would contribute to the axial Ward identity. Up to ^{g^) these 
parity violating terms arising out of {x"^(,x)) are as follows, 

(/) : ~2im Trace W^"' 

vanishes because of Dirac Trace 

-2im Trace W^^'Vi^^'^Yi^'" - 2im Trace W';f'"Vf'if'"Vf^'" . 

(12) 



Contribution to the anomaly from (x'^ix)) can be calculated by making the follow- 
ing changes from the standard Wilson case: 



and 



Trace ys{W +m)'^Vi'^Vi<:^W Trace y5WT"Vi<^'"Vi^"'W 



Trace Y50^Vi'i^[0,W]'^W Trace 750^'"yi^'" [0, 



(13) 



(14) 



Thus compared to the expression in the case of standard Wilson term, the expression 
for anomaly in the case of Osterwalder-Seiler Wilson term has the following features: 
(i) absence of fermion mass m in the numerator and (ii) absence of mixing between 
mass term and Wilson term (Win) in the denominator Explicitly, in place of Eq. (|6]l 
we find the contribution to FF from (x'^ix)), in the case of the Osterwalder-Seiler 
twisted Wilson term 

{X"') = e^ypiF^v{x)Fpx{x)Y,cos{p^a)cos{pva)cos{ppa) 



X Wo{p) cos{pxa)Woip)- 4rsm^ip^a) {%'{p)f 



(15) 



In the following we refer to the terms proportional to cos{pia) and sin^(p;^a) as DD 
and DW terms respectively since they correspond to contributions from ViVi and ViVj"* 
terms in the trace. 

In the Axial Ward Identity, the mass term can be written as 

-2ot Trace = -2m Trace 75(0- m - /75W)(^'" 

+^''W"y^'" + •••)■ (16) 

'•^), the parity violating contributions from the mass term are 

+2im Trace W^"' 
-2im Trace W]^'"' 

+2im Trace W^""yi^'"Yi^'"' + 2/m Trace . 

(17) 
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Figure 1: The function I{am,r,L)foi r = I for unimproved, ^(fif) improved, G{c?-) improved, and OStm 
Wilson fermions for tlie range of am between 0.01 and 1.0 at L = 40. 



Comparing Eqs. ( 12 1 and ( 17 1, we immediately see that and ^{g^) parity vio- 

lating terms cancel between mass term and 
The ^{g^) term in Eq. ([Tt} 



ism{ap^)){cos{app) 
i sm{a p^)) {cos {a pp) 



sm{ap^)Ffip (cos(fl/?^) 

-(cos(fl/?^) 
— (cos (fl ) + / sin(a/)^ ) ) (cos {a pp ] 

+{cos{apf^) — i sm{a Pn)) (cos {a Pp) 



-ism{app)) 

- ism{app)) 

- ism{app)) 

- ism{app 



0, on summation. 



(18) 



Thus we explicitly verify, in the case of Osterwalder-Seiler twisted Wilson term, the 
cancellation of parity violating terms in the flavour singlet Axial Ward Identity up to 

^{g^l 
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3. Comparison with unimproved, ff{a) and ff{a^) improved Wilson fermions 

Consider the following next to nearest neighbour interaction term added to modify 
the standard Wilson term JSl: 



AS' 



\lf{x)Ux,^Ujc+fi.jiW{x + 2n) 



Xj/ix + 2ii)ul+^^^Ul^\li{x) - 2^t{x)^i{x) 



(19) 



The coefficient ^ is chosen so as to cancel ff{d) contributions to the tree level fermion 
propagator and the fermion-gluon vertex. Now the total fermion action is 5 = 5/r + 
AS^ . From the corresponding flavour singlet axial Ward identity, we find the total 
contribution to the axial anomaly 



2£)ivpX Fiiv{x)Fpx{x) Y,cos{ap^)cosiapp)cos{apx) [^o(/?)]^ 
p 

cos{apy)[m + Wo{p) + W^ip)] 

-4rsin(flpv)(sin(flpv)-^sin(2flpv))][Wo(p)+Wo^(/7)] . (20) 



Here 



Woip)+Wiip) = I £ y-il-cos{ap^)) + ^(-1 +cos{2apf,)) 



f ^Lsi"^(«/'M) + ("^ + L-[l"'^°s(«/'M)] + J-[-l+cos(2fl/7^)])M . (21) 



Next consider the modification of the kinetic term with the next to nearest neigh- 
bour interaction ll8ll9l[T0l: 



1 

16a 



(22) 



so that the total fermion action S = Sf + AS' + AS" is 6{c?-) improved and leads to the 
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total contribution to the chiral anomaly 

{X") = ^£^vpX i F^y{x)Fpx{x) ^nl^i(cos(fl/,;L) - ^cos(2flM)) K{P)? 

3 



3--(Wo(p)+Wo^(p)) 



-m + 3 

4 a a' 



{cos{ap^) - \ cos{2ap^))^ + {-ia{D + D')o)n{p){sin{ap^) - i sin(2a/?^)) 
(cos(a/7^)- icos(2a/7^)) 



with {D + D')ox{p) = j^[sm{pxa) - lsm{2pxa)] and 



(23) 



^Um + D')o)^){p)}^ + {lm + Wo{p)+W^{p)y 



-1 



(24) 



Fig. [T] compares the function I{am,r,L) for r = 1.0 for unimproved, &{a) im- 
proved, ff{a^) improved and OStm Wilson fermions for the range of am between 0.01 
and 1 .0 at L = 40. Note that for clarity, am is plotted in the logarithmic scale. The 
figure shows that in the range of am studied, G{cp-) corrections are relatively small 
compared to the &(a) corrections and the approach to the chkal limit is the fastest and 
the flattest for the OStm Wilson fermions. 



4. Numerical Results for Osterwalder-Seiler twisted mass Wilson Fermion 

In this section we present numerical results for the function /(am, r. L) for Osterwalder- 
Seiler twisted mass Wilson fermions. Since our main concern is the approach to the 
infinite volume chiral limit in the continuum, we study in detail the quark mass (am), 
the Wilson parameter (r) and the finite volume (L) dependence. Since a positive semi- 
definite transfer matrix is guaranteed only for r < 1, we restrict our study to the range 
< r < 1. 



4.1. Quark mass dependence 

In Fig. |2]we show the function I {am, r, L) for OStm Wilson fermion for the range of 
am between 0. 1 and 1 .0 at L = 40. We find that for most values of am, DW term over- 
shoots the answer, which is compensated by the DD term. Cutoff effects are present in 
the mass range am > 0.1, but the result is very close to the continuum chiral answer for 
am < 0.1. 



4.2. Wilson parameter r dependence 

For unimproved, ff{a) and &(c^) improved Wilson fermions, it has been demon- 
strated that in the chiral limit, the function I{am,r,L) is independent of the Wilson 
parameter r. For the OStm Wilson fermions, it is also interesting to check the r depen- 
dence for non-zero am. In Fig. |3]we show the r dependence of terms with integrands 
proportional to cos{pia) (DD) and sin^{pxa) (DW) in Eq.( 15 1 and the sum of the two 



7 



1 • 



- 



# DD Term 
■ DW Terra 

• Total 



4. ^^♦♦^ 



0.1 

am 



Figure 2: The function I{am, r,L) for r = 1 for OStm Wilson fermions for tlie range of am between 0.01 and 
1.0forL=40. 



contributions to the function /(am, r,L) (DD and DW) for am = .01 (right) and am = .1 
(left) for OStm Wilson fermions at L = 80. We note that the first term (DD) is dominant 
at small r and the second term (DW) is dominant at large r. Even though the two indi- 
vidual contributions have very strong r dependence, the sum is seen to be independent 
of r to a good numerical accuracy for am = .01. 

4.3. Finite volume dependence 

In Fig. [4] we compare the finite volume dependence of the function I{am,r,L) for 
OStm Wilson fermions at r = 1 and range of am between 0.01 and 0.1. In this range of 
am, convergence is satisfactory as volume increases but the convergence rate becomes 
slower if the value of am is smaller. 

In conclusion, we have shown that parity violating terms do not arise in the flavour 
singlet axial vector Ward identity up to ff{g^) for the OStm Wilson fermions and the 
approach to the chiral limit in comparison with (^{a), and Gic?-) improved Wilson 
fermions is very satisfactory. 
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0.01 (right) and 0.1 (left) at L — 80. The contributions from DD and DW terms are also shown separately. 
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Figure 4: The function /(am, r,L) for OStm Wilson fermions for the ranges of am between 0.01 and 0.1 for 
different L. 
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